Force networks and the dynamic approach to jamming in sheared granular media 
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Diverging correlation lengths on either side of the jamming transition are used to formulate 
a rheological model of granular shear flow, based on the propagation of stress through force chain 
networks. The model predicts three distinct flow regimes, characterized by the shear rate dependence 
of the stress tensor, that have been observed in both simulations and experiments. The boundaries 
separating the flow regimes are quantitatively determined and testable. In the limit of jammed 
granular solids, the model predicts the observed anomalous scaling of the shear modulus and a new 
relation for the shear strain at yield. 



Jamming takes place when an amorphous collection 
of particles spontaneously develops rigidity and supports 
weight like a solid instead of flowing like a liquid 
The transition occurs without static spatial ordering, 
but is accompanied by long range dynamical correla- 
tions arising from the collective motion of groups of par- 
ticles 0, S Hi- I n the case of sheared granular materials, 
the transition from jammed to flowing phases occurs as 
the applied stress S is increased or the packing fraction 
4> is decreased, and there is a well defined packing frac- 
tion 4 at which the system jams in the limit of zero 
stress [H, 0| . Diverging correlation lengths are observed 
on each side of the transition, and are related to the av- 
erage size of force chain networks for <p < <fi c 0, 0, H and 
the average size of isostatic clusters for tf> > 4> c [3|, |9| . 

In the jammed state many macroscopic observables ex- 
hibit power law scalings in {4>— <p c ) [H, fi[ and the flowing 
state rheology changes dramatically as the packing frac- 
tion is increased above cj> c . Although the jamming tran- 
sition controls both dynamic and static properties near 
<f) c , theories tend to focus mainly on the latter 0, [13, 12 1. 
In this Letter we theoretically explore the jamming tran- 
sition by first studying flows with non-zero shear rate 7 
and then taking the limit of 7 — > to access the static 
case. This procedure leads to quantitative predictions for 
the flowing rheology that match observations in Ref . [II] • 
It also predicts the anomalous static scaling of the shear 
modulus measured in Ref. [5, 6] and a new scaling rela- 
tion for the static yield strain, 7* oc {<j> — 0c) 1 , which 
provides a testable prediction of the theory. 

Granular materials behave as peculiar liquids. This 
can be clearly demonstrated by measuring the stress ten- 
sor S. Depending on various parameters, the system will 
exhibit either Bagnold scaling S a j 2 , elastic-inertial 
scaling £ cx 7 1 , or quasi-static scaling S oc 7 (i.e. con- 
stant) [TH . The rheology of the flow is sensitively depen- 
dent on 4>, with Bagnold scaling for <f> < <p c , quasi-static 
scaling for cf> > <p c , and elastic-inertial scaling in between. 
This is in marked contrast to Newtonian fluids where 
S cx 7 and only the proportionality constant depends on 
(j). Although the phase diagram of granular shear flow has 
been extensively studied in simulations [3 H El and 
experiments |15[ , the origins of the rheological crossovers 



remain to be explained. Arriving at a solution to this 
problem requires a well developed theory for the stress 
tensor in dense granular flows. 

The stress tensor ultimately depends on the amorphous 
microscopic arrangement of grains. Forces are transmit- 
ted via contacts between grains and a perturbation on 
one grain can have long range effects. Indeed, both ex- 
periments and simulations indicate that contact forces 
are correlated 0, EH and tend to form quasi one dimen- 
sional filaments, or force chains, that permeate the mate- 
rial These force chains assemble into network struc- 
tures with a characteristic size £ that diverges as <f> — > <f> c 
from below @, 0, @|, as illustrated in Fig. [TJ Since con- 
tact between grains is the only form of force transfer, S is 
fully determined by properties of the networks @, El, • 
Here we investigate the role of force chain networks in de- 
termining the value of the stress tensor by constructing 
a model of momentum transfer. 

The Force Network Model (FNM): The central concept 
of the FNM is that the force F i3 between a pair of con- 
tacting grains {i,j} can be expressed as the sum of a 
colhsional part F? 3 and a static part F* 3 ' . The collisional 
part is proportional to the square of the relative velocity 
between grains i and j upon initial contact, is well de- 
scribed by kinetic theory [2(J, and scales with 7 2 . It is 
the force that would be expected from collisions between 
pairs of grains. The static force F* 3 is a consequence of 
the pressure induced by the network surrounding a pair 
of grains Collisional forces from other contacts in 

the network are transferred to from first nearest 

neighbors, second nearest neighbors, and so forth, all the 
way to the edge of the network. The presence of multi- 
ple contacts in the cluster thereby increases the force at 
and this contribution can be expressed as a sum 
over all paths between the contact and every other 

contact {i',j'} in the connected force network: 



f: 3 



NiQiFl 
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In this equation, £ is the path-length (in grain diame- 
ters) between the two contacts {i, j} and {i',j'}, Ng gives 
the number of pairs of contacts connected over a path of 
length £, and F^ 3 is the collisional force on the contact 
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at the end of the path. The function Qg accounts for the 
fact that only a portion of the collisional force is trans- 
ferred at each link in the path, and no more than the 
collisional force can be transferred from each contact in 
the network. In constructing this equation we have as- 
sumed that grains are perfectly rigid, which ensures that 
the propagation of collisional forces is instantaneous. 

The FNM extends previous theories of granular flow 
that separate the stress tensor into static and collisional 
components plj by introducing a natural physical ba- 
sis for calculating the static stress. Static stresses arise 
from clusters (not necessarily percolating) of simultane- 
ously contacting grains, which exist even in the limit of 
infinitely rigid grains @, 0, Q • These force networks are 
created by dissipative contacts between grains and de- 
stroyed by the shear flow. Since each of these rates is 
proportional to 7 in the rigid grain limit, the properties 
of the steady state force networks are independent of 7. 

Bagnold Scaling: Eqn. fl} yields a constitutive relation 
for each independent component of the stress tensor that 
does not contain any adjustable parameters and is ex- 
pected to hold over scales much larger than the network 
size £. This prediction has been shown to match data 
from simulations for all packing fractions <p < 4>c @> Ell • 
A major feature of the constitutive relations is that Bag- 
nold scaling holds, or S ocy 2 . This is because all of the 
network parameters in Eqn. ([1]) are independent of 7 and 
the rheology is thereby set by the average collisional force 
between grains, as predicted by kinetic theory (20j . 

The form of the constitutive relations is especially sim- 
ple for rigid granular materials slightly below (f> c . At large 
packing fractions, each grain has many contacts and the 
entire collisional force from every contact is transferred to 
its nearest neighbors. Thus Q\ equals its maximum value 
1/Ni, which ensures that Qi = 1/Ng for all values of I. In 
this limit, it follows from Eqn. fl} that S ~ £(-Ff> c ) oc £7 2 . 
The stress tensor is predicted to depend linearly on the 
size of the force networks, independent of spatial dimen- 
sion, and is insensitive to other network parameters. 

The Elastic- Inertial Regime: The FNM prediction for 
the stress tensor has been constructed for perfectly rigid 
grains, but must be altered for realistic grains with a fi- 
nite stiffness. In particular, excitations do not propagate 
instantaneously through force networks, but occur at a 
finite speed c, which is an increasing function of the stiff- 
ness of the grains. Combined with the average lifetime 
of the force chain networks r, the speed of force propa- 
gation defines a length scale cr that gives the maximum 
extent forces can propagate through networks. 

If £ < cr then forces will propagate through the en- 
tire force network before it is destroyed. Therefore the 
assumptions of the rigid grain FNM hold and X oc £7 2 , 
averaged over time scales larger than £/c. However, if 
£ > cr, the collisional force from a single contact will not 
be transferred to the entire network, but only to contacts 
a distance cr away. Thus the network size £ is no longer 
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FIG. 1: A schematic illustration of three important length 
scales, plotted as a function of packing fraction <j>, and their 
relation to granular rheology near <f> c . £ gives the size of force 
chain networks and diverges as (f> — > I* quantifies the size 
of marginally stable regions in the material and diverges as 
4> — > (f>c . Over length scales larger than Z*, the system re- 
sponds as a conventional solid, cr represents the maximum 
correlation length through force networks and increases with 
packing fraction. The macroscopic rheology (indicated be- 
neath the graph) depends on the smallest length scale, and 
transitions occur at the intersections. <f) c is the packing frac- 
tion where the system jams in the limit of zero stress. Our 
predictions for 0h s and qs are given in Eqns. (J5j [6}. The 
numerical values of the length scales depend on many param- 
eters. Here we have used £ = (<f> c — <f>)~ 3 ^ 2 , I* = (<j> — 4> c )~ 2 
and cr = 10 + 3O0 3 for visualization, but the crucial features 
are the intersections rather than the specific numerical values. 

relevant and the size cr of the correlated networks should 
be used instead. Since £ diverges at <f> c , there will always 
be a critical packing fraction </>h s above which cr < £, as 
illustrated in Fig. [TJ For <j> > tfihs the shear flow will be 
in a regime where £ must be replaced by cr. 

The value of c is a monotonically increasing function 
of the stiffness of grains and is also weakly dependent 
on the geometry of contact networks. It is a property of 
the granular material and does not depend on the shear 
rate. The average network lifetime r is proportional to 
7 _1 since the shear flow pulls apart the networks 12, 13j |. 
Writing r = 77/7 and substituting £ — > cr yields X oc crjj. 
This predicts a linear scaling with 7, as was observed in 
Ref. [l2j and named the elastic-inertial regime. 

The introduction of finite grain stiffness to the FNM 
produces a crossover from Bagnold's scaling to elastic- 
inertial scaling as the packing fraction or grain stiffness 
is increased. The elastic-inertial scaling reveals that a 
second time scale, shown here to be proportional to the 
speed of force propagation through force networks, be- 
comes important at high packing fraction. 

The onset of quasi-static flow: For granular materials 
with <f> > <p c , another length scale I* becomes relevant Q 
that is related to the departure of the material from the 
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isostatic limit. Isostatically jammed materials are con- 
figured such that contact forces can be determined iden- 
tically from the constraint that no particle moves [22I ]. 
This occurs when the coordination number z, equal to 
the average number of contacts per particle, approaches 
a critical value z c that depends on the spatial dimension. 
Contact forces are highly correlated in the isostatic state 
since breaking any contact will produce a cascade that 
alters the value of every other contact force. 

The length scale I* quantifies the extent of cascades 
that result from contacts being broken. For a system 
with coordination z = z c + Sz, a disturbance that breaks 
a single contact will result in a cascade over a length 
scale of I* oc l/fo @, a scaling that has been verified 
in simulations [3J]. The material is therefore fragile over 
length scales smaller than I*, but the excess contacts 5z 
stabilize disturbances over length scales greater than I*. 
The value of I* is also related to the packing fraction, 
since Sz oc (</> - ^c) 1 / 2 di, and is plotted in Fig. [TJ 

The presence of the length scale I* has important dy- 
namic effects when combined with the FNM and the 
length scale cr. When cr < I*, forces propagate over 
regions that are effectively isostatic and only marginally 
connected. In this limit, if a contact is broken then there 
are rearrangements of grains over the entire range of the 
force chain networks and inertial scalings are important. 
However, if ct > I*, forces propagate over distances that 
are large compared to the rearranging regions, redundant 
contacts stabilize the flow, and forces are no longer in- 
ertial. Thus, when ct = I* the system crosses over to a 
quasi-static regime with £ oc c 2 /I* for ct > I*. 

Predictions at non-zero shear rate: The above consid- 
erations characterize the dependence of the stress tensor 
on the length scales £, ct, and I*. This prediction is ob- 
tained using the FNM for rigid grains and carried through 
to the elastic-inertial and quasi-static regimes by investi- 
gating crossovers in the characteristic length scales. The 
stress tensor is given by 

C £,j 2 , for £ < ct; 

£ oc < cr/j, for £ > ct and I* > ct; (2) 
I c 2 r/ 2 /l*, for £ > ct and 1* < ct; 

and a schematic of the various length scales and flow 
regimes is shown in Fig. [TJ 

Instead of characterizing these regimes in terms of 
length scales, as is done in Eqn. {2]), it is advantageous 
to recast the transition points in terms of a critical shear 
rate. The transition from Bagnold to elastic-inertial scal- 
ing occurs when £ = ct. Given that t = 17/7, we can 
define a critical shear rate A f^ s . Similarly, the transition 
from elastic-inertial scaling to quasi-static scaling occurs 
at a critical shear rate 7* s . The critical values are 

7hs = c?7/£, (3) 

and 

%* s = cn/r. (4) 
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FIG. 2: The predicted flow-map of granular materials, as 
a function of packing fraction <f> and normalized shear rate 
7' = 7/077. The stress tensor scales with -y 2 in the Bagnold 
regime, with 7 1 in the elastic-inertial regime, and with 7 in 
the quasi-static regime. The scaling behavior of <f> qs and (f>hs 
are taken from Eqns. (|5I6[1 . using ip = 1. 

Major features of the resulting granular rheology are il- 
lustrated in Fig. O 

For a granular shear flow with constant <fi < <j> c , cor- 
responding to a horizontal slice through Fig. [21 the sys- 
tem exhibits Bagnold scaling for 7 < 7^ s and elastic- 
inertial scaling for 7 > 7^,. This unexpected behav- 
ior where Bagnold's scaling, normally associated with 
"rapid" flows, actually occurs for small 7 in flows with 
constant packing has been observed previously fl2| . For 
constant 4> > <p c , the system exhibits quasi-static scaling 
for 7 < 7* s and elastic-inertial scaling for 7 > 7* s . The 
emergence of quasi-static flow as the shear rate is reduced 
in dense materials is a feature that has been observed in 
experiments [l5| and simulations [23| . 

The flow regimes can also be characterized by the pack- 
ing fractions </>h s and 4> qs at which the transitions occur. 
For ct > 1, Fig. [JJ guarantees that the difference be- 
tween the critical packing fractions and <p c will be small. 
We can therefore use the scaling form I* oc (cj) — 4 > cY~ 1 ^ 2 
for the isostatic length scale in the dense regime 0, 0] . 
Combined with Eqn. this predicts 

4> qs - <j) c oc (— ) 2 . (5) 
cq 

Thus as 7/ cq approaches zero, the crossover to quasi- 
static flow occurs closer to <p c . The functional form of 
4> qs is plotted in the flow-map of Fig. [2] 

We can also predict the value of 0hs in a similar way. 
Since £ diverges at c> c , then it must scale as £ oc (</> c — 4>)^ 
near the transition. For ct> 1 Eqn. ^ predicts 

t 7 \ 

<j) c - (f> hs oc (-j . (6) 

Eqns. (|5l6p fully determine the crossovers in rheology and 
allow us to construct the flow-map in Fig. O This figure 
is consistent with the numerically obtained flow-maps of 
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granular rheology in Fig. 10 of Ref. using c oc 
where k is the grain stiffness. 

Predictions at zero shear rate: In addition to predicting 
the rheology of granular materials, Eqn. ([2|) also gives 
the scaling of the stress tensor at zero shear rate when 
the system is jammed. In the limit of 7 — * the yield 
pressure p y and shear stress s v are both proportional to 
c 2 /l*. For stresses less than these values, the system [4 
behaves elastically with s = G7, where G is the shear 
modulus and 7 is the shear strain. We can therefore set 



Vk, 



s y = Gj* , with 7* the critical shear strain at which the 
system yields. The speed of force propagation c is then 
given by the speed of elastic shear waves: c cx \[G. 

The dependence of G and p y on packing fraction (or 
excess coordination) in the jammed state has been the 
focus of much attention. The behavior of p y can be un- 
derstood by assuming affine grain motion, whereas the 
scaling of G is anomalous and affected by particle re- 



arrangements after a strain step 0. 
Eqn. |2]) predicts that 

G oc py/Sz, 
7* oc Sz. 



Using I* 



and 



(7) 
(8) 



The first scaling relates the shear modulus to the yield 
pressure and excess coordination. It matches measure- 
ments fromsimulation for all dimensions and interaction 
potentials [!, [|| . The second relation predicts the critical 
shear strain at which the system yields. Although this 
has not yet been measured in simulations or experiments, 
the FNM predicts that yielding is a purely geometric phe- 
nomena, unaffected by interaction potential. 

Conclusions: We have taken a dynamic approach to 
the jamming transition by using the FNM to predict the 
stress tensor in the flowing regime and then interpolat- 
ing to the jammed state. While this procedure was aimed 
at understanding athermal granular materials, a similar 
technique should apply in amorphous systems that jam 
as the temperature is reduced. Since the isostatic length 
scale £* exists for all jammed amorphous materials, the 
challenge is to identify relevant correlation lengths in the 
flowing regime and relate them to properties of inter- 
est. Jamming may then be universally understood as the 
transference of correlation from flowing length scales to 
the isostatic length scale. 
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